Andreev drag effect via magnetic quasiparticle focusing in normal-superconductor 
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We study a new hybrid normal-superconductor (NS) 7r-junction in which the non-local current 
can be orders of magnitude larger than that in earlier proposed systems. We calculate the electronic 
transport of this NS hybrid when an external magnetic field is applied. It is shown that the non-local 
current exhibits oscillations as a function of the magnetic field, making the effect tunable with the 
field. The underlying classical dynamics is qualitatively discussed. 
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Electron-transport properties of normal- 
superconductor hybrid nanostructures have been 
the subject of extensive theoretical and experimcntali 
attention. Experiments carried out on nanostrucures 
containing fcrromagnets (F) and superconductors 
(S) reveal novel features, not present in normal- 
metal/superconductor (N/S) junctions, due to the 
suppression of electron-hole correlations in the ferro- 
magnet. When spin-flip processes are absent, further 
effects are predicted, including the suppression of con- 
ventional giant magnetoresistancc in diffusive magnetic 
multilayers^ and the appearance of non-local currents 
when two fully-polarized ferromagnetic wires with 
opposite polarizations make contact with a spin-singlct 
superconductor—. The latter effect, also called the 
Andreev drag effect, has been highlighted, because of 
interest in the possibility of generating entangled pairs 
of electrons at an N-S interface^. A recent study of 
such a junction in the tunneling limits^ predicts that the 
magnitude of the non-local current decreases exponen- 
tially as exp (— 2L/7r£ c ), where £ c is the superconducting 
coherence length, and L is the distance between the F 
contacts. The effect can be enhanced by inserting a 
diffusive normal conductor between the superconductor 
and the ferromagnetic contacts leads as shown in Ref&i. 
The off-diagonal conductance^, which is always negative 
in the normal case, can have a positive value of order 
the contact conductances of these systems. However, 
the value of the off-diagonal conductance is determined 
by fixed material parameters, such as the polarization of 
the ferromagnets, and the spin-flip time in the normal 
diffusive conductor. Therefore, it is of interest to study 
alternative methods for material-independent tuning of 
the non-local current. 

In this work, we show that even in the absence of fer- 
romagnetic contacts, an enhanced Andreev drag effect 
is possible with the N/S structure shown in Fig. ^ We 
shall demonstrate that the non-local current is enhanced 
by orders of magnitude compared with the structure in 
which ferromagnetic leads were used to detect the cur- 



rent—. Moreover, the magnitude of the non-local current 
can be tuned by varying a magnetic field applied perpen- 
dicular to the system. The necessary field is much lower 
than the critical field of the superconductor. 
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FIG. I: (color online) The hybrid N/S nanostucture con- 
sists of an infinitely long N/S ballistic waveguide, comprising 
a normal (N) metal of width Wn coupled to a spin-singlet su- 
perconductor (S) region of a width that is much larger than 
the superconducting coherence length £ c . To measure the 
non-local current, two ballistic normal leads, at voltages vi 
and V2, and with width Wl, separated by a distance L, are 
in contact with the normal waveguide. The left and right 
ends of the waveguide act as drains which absorb any quasi- 
particles exiting to the left or right. The magnetic field B is 
applied perpendicular to the system (in our calculation B > 
corresponds to a field pointing out of the plane of the system). 

To calculate the non-local current we employ 
the current-voltage relation developed for nor- 
mal/superconducting hybrid structures in the linear 
response limit—. Assuming that the voltages vi at lead 2, 
is the same as the voltage v of the condensate potential, 
for the arrangement shown in Fig. ^ one finds that the 
currents in lead 1 and 2 are 

2e 2 

h = — (N-R + R a ) (vi - v), (la) 



2e 2 

h = — (T a - T ) («! - v), 



(lb) 



where v\ is the voltage at lead 1 and is the number 
of open scattering channels in the normal leads of width 
Wl ■ Here i?o (To ) are the reflection (transmission) coef- 
ficients for an electron from lead 1 to be reflected (trans- 
mitted) to lead 1 (2), and R a (T a ) are the Andreev 
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reflection (transmission) coefficients for an electron from 
lead 1 to be reflected (transmitted) to lead 1 (2) as a 
hole. R a and Ro satisfy the inequality N — R + R a > 0, 
thus, I\ is always positive for positive v\ — v. All coeffi- 
cients are evaluated at the Fermi energy using an exact 
scattering matrix formalism. 

It is easy to see from Eq. Q that whenever Andreev 
transmission dominates normal transmission (ie T a > To) 
the currents I\ and I 2 have the same signs, ie a current 
in lead 1 induces a current in lead 2 flowing in the same 
direction. In semi-classical point of view, this means that 
hole like quasi-particles leave the system at lead 2. This 
is the Andreev drag effect. On the other hand, in the case 
when the normal transmission is larger than the Andreev 
transmission (ie T a < To), electron-like quasi-particles 
leave the system through lead 2, yielding a current flow- 
ing opposite to the direction of the current in lead 1. 

In what follows now, we show that for the system de- 
picted in Fig. [2 the ratio T a /T can be tuned by an ap- 
plied magnetic field. To this end, we calculate the trans- 
mission coefficients for the system using the Green's func- 
tion technique^ developed for discrete lattice. Each site 
is labelled by discrete lattice coordinates (x, y) and pos- 
sesses particle (hole) degrees of freedom ip e ( h \x, y). The 
magnetic field is incorporated via a Peierls substitution. 
In the presence of local s-wave pairing described by a su- 
perconducting order parameter A(x,y), the Bogoliubov- 
de Gennes equationiS (BdG) for the retarded Green's 
function takes the form 
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where the components of H arc 
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(2b) 



Here Ep is the Fermi-energy, and n x and n y are the near- 
est neighbors of (x, y) in the x and y direction, respec- 
tively. 

Within the Landau-gauge with a vector potential in 
the x-direction, ^ x = 7oe* 9 ^, ^ y = 70, where 70 is the 
hopping parameter without magnetic field. The phase 
9{y) for the Peierls substitution is zero in the supercon- 
ducting region, and it is given by 9(y) = Ba 2 7r(W?f — 
y)/&o m the normal region, where a is the lattice con- 
stant™, and $0 = h/2e is the flux quantum. This choice 
of gauge results in a uniform magnetic field B in the 
normal region, and zero magnetic field in the S region, 
while the translation invariance in the x direction is pre- 
served. The order parameter is assumed to be a step 
functior*i2ii£, ie constant Ao in the S region and zero 
otherwise. The phase 9 is set to #i oa d = Bci 2 itWn /&o 
in the leads 1 and 2 to ensure the continuity of the vec- 
tor potential. The parameters of the Hamiltonian H are 



chosen to model an experimentally-realizable situation in 
the quasiclassical regime, ie Wn ^> Fermi wavelengthAi. 

From the Green's function and the scattering matrix 
for the system, the transmission and reflection coeffi- 
cients are calculated as a function of the magnetic field. 
Our central result, shown in Fig. |3 is that the differ- 
ence between the Andreev and normal transmission co- 
efficients T a — To (which proportional to the measurable 
current I2 according to Eq. Ijlb|) ') is an oscillating function 
of the magnetic field. Furthermore, since positive peaks 




FIG. 2: (color online) T a — To as a function of the magnetic 
field (in units of $o/(2a 2 7r)) at the Fermi energy Ef- In lead 
1 only one mode was allowed. The wave functions at magnetic 
fields corresponding to letters A and B on the peaks of the 
curves will be shown in Fig. 

correspond to pronounced Andreev drag effect and the 
heights of the positive peaks are comparable with those 
of the negative peaks, the non-local current can be as 
large as the normal current in our hybrid system. 

A striking feature of Fig. |3 is that it is an asymmetric 
function of B. This can be understood qualitatively by 
tracing the classical cyclotron orbits of quasi-particles, 
bearing in mind that when electron-hole conversion oc- 
curs at the NS boundary, the chirality of the electron- 
like and the hole-like orbits is preserved and therefore a 
geometrical construction for their classical trajectories is 
different from that of normal systems^. Examples of tra- 
jectories obtained from this new geometrical construction 
are plotted in Fig. El 




FIG. 3: (color online) Classical trajectories in a magnetic 
field. Figure (a) shows the case when the electron does not 
reach the N-S interface. New types of trajectories involving 
Andreev reflections are sketched on figures (b) and (c), for the 
cases when only the electron (b), both the electron and the 
hole (c) can reach the side of the waveguide where the leads 
are attached. Blue (red) lines refer to the electron (hole). 
Electrons are injected perpendicularly into the waveguide at 
the positions marked by arrows pointing up. 
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For B > electrons injected from lead 1 will follow 
classical orbits bending to the left. For a large enough 
\B\ these will exit to x — — oo, without reaching the N- 
S interface, and impinging on lead 2, as can be seen in 
Fig. Therefore for large positive B all transmission 
coefficients from lead 1 to lead 2 vanish. Andreev re- 
flection can occur if \B\ is sufficiently small to allow the 
electrons to reach the N-S interface. This condition is 
defined by \B\ < Bi, where B\ is the field for which the 
cyclotron radius R c = Wn, where R c = y/2mEp/eB. As 
shown in Figs. |3J>c, the transport direction is reversed 
compared to the normal case due to Andreev scatter- 
ing, because even if the classical orbits are anti-clockwise, 
quasi-particle transport is to the right, resulting in quasi- 
particles impinging on lead 2. This is why the asymmetry 
in Fig. |21 arises. On the other hand, as shown in Fig. 0t>, 
if R c is not sufficiently large, there is no drag effect, be- 
cause the trajectories of the holes do not hit the side of 
the waveguide to which the leads are attached. 

Andreev drag effect only occurs for \B\ < -B max , where 
the maximum field -B max is determined from the condi- 
tion R c > 2 Wat. By appropriate choice of the width Wn 
of the normal part of the waveguide, -B max can be much 
less than the critical field of the superconductor. The 
trajectory relevant for this case is shown in Fig. On 
the normal side of the waveguide, normal quasi-particle 
reflections alternate between electrons and holes, sepa- 
rated by equal distances 5. Assuming that the electrons 
are injected perpendicularly into the waveguide, simple 
geometrical considerations give the following condition 
for maxima in T a : 

L = (2n + 1)5, (3) 

where n is an integer counting the number of normal re- 
flections of the hole at the side of the norm al wavegui de 
to which the leads are attached, and S = 2^J R? c — W N — 
y/RZ - W%-R c . FromEq. © one can find a magnetic 
field B n for each n. The peaks in T a can be expected at 
B n . Taking into account the finite widths of the two leads 
we calculated the ranges of B for each n in which a peak 
in T a should be found, which corresponds to the range 
of B for which a classical trajectory of the hole hits the 
finite-width interface of lead 2. In Fig. 01 we plotted the 
ranges of B n as vertical bars together with those values 
of magnetic field at which we obtained peaks from the 
exact quantum calculations. One can see that the agree- 
ment between the quantum and the classical calculation 
is reasonable. Thus, this simple classical argument can 
be used to estimate the magnetic field needed to obtain 
enhanced Andreev drag. 

To reinforce the above detailed classical picture, we 
now calculate the electron and hole component of the 
wave functions inside the N/S waveguide, and compare 
them with classical orbits. The contribution of the n-th 
incoming mode (from the left N-lead of width Wl) to the 
wave function at point rg = (x, y) of the waveguide is 

r n W (r s ) =J2G e < he \r s ,v L ) X +' e (r L ), (4) 
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FIG. 4: The range of B n shown as a vertical bar for each 
n, together with those values of magnetic field at which we 
obtained peaks in T a from the quantum transport calculation. 
(Green line is connecting these peaks, but only for guiding the 
eyes). 

where runs over the surface of the left lead. Here 
the appropriate components of the retarded Green's func- 
tion are defined in Eq. Q2a|). and xi' e { r L) is the trans- 
verse wavefunction of the n-th incoming electron channel 
of the left lead, normalised to unit flux. The modulus 
square of the electron and hole components of the wave 
function are shown in Fig. [5] for two different magnetic 
fields, corresponding to the positive peaks A and B in 
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FIG. 5: (color online) From top to bottom the electron (left) 
and hole (right) probability amplitudes are plotted, corre- 
sponding to the Andreev-transmission peaks marked by A and 
B in Fig. respectively. In our geometry Wn = 80, W = 100, 
lead 1 (lead 2) is located at -10 < x < (100 < x < 110). 
Distances are in units of the lattice constanii— . 

Fig. |21 For these scattering states, the hole probability 
amplitude has a local maximum at lead 2. There are sev- 
eral other maxima of the probability amplitudes of the 
wave function at the lower side of the waveguide, both 
for the holes and for the electrons. For each positive peak 
of (T a — T ), the condition Eq. (J3J is satisfied, where n 
is the number of maxima of the hole probability ampli- 
tude between the leads, and 5 is the distance between the 
nearest electron and hole maxima. 
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In conclusion, we have shown that even in the absence 
of ferromagnetic leads, an enhanced non-local current can 
be obtained by including a normal region between the 
leads and superconductor, and applying magnetic fields 
perpendicular to the system. The current flowing from 
lead 1 to lead 2 shows oscillations with alternating signs 
as a function of magnetic field in the small-field regime, 
corresponding to alternating magnetic focusing of elec- 
tron and hole-like quasi-particles between the two leads. 
Unlike an earlier proposal^, where T a is exponentially 
suppressed with lead separation, the non-local current 
remains significant even for a lead separation much big- 
ger than the coherence length of the superconductor. We 
discussed how the quantum results could be interpreted 



qualitatively in a fully classical treatment providing a 
better insight into the Andreev drag effect in our system. 
For the future it would be of interest to extend the semi- 
classical approach developed for normal focusing prob- 
lem. In this analysis one has to involve the semi-classical 
wave functions of the particles taking into account the 
more complex caustics^ formed for both electrons and 
holes. 
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